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We provide a model-independent determination of the quantity Bo(md — m u )- Our approach rests 
only on chiral symmetry and data from the decay of the eta into three neutral pions. Since the 
low-energy prediction at next-to-leading order fails to reproduce the experimental results, we keep 
the strong interaction correction as an unknown parameter. As a first step, we relate this parameter 
to the quark mass difference using data from the Dalitz plot. A similar relation is obtained using 
data from the decay width. Combining both relations we obtain Bo(md — m u ) = (4495 ±440) MeV 2 . 
The preceding value, combined with lattice determinations, leads to the values m u (2GeV) = (2.9 ± 
0.8) MeV and m d (2 GeV) = (4.7 ± 0.8) MeV. 

PACS numbers: 12.39.Fe, 11.30.Rd, 14.40.Aq 
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INTRODUCTION 



One of the most reliable ways to extract predictions in 
the low-energy regime of strong interactions is the use of 
the chiral symmetry of QCD. This technique, suggested 
in [l[ , was put on ground in 0, H[ and has been fully devel- 
oped in chiral perturbation theory over the last 30 years. 
In the pion sector one finds strong evidence for the cor- 
rectness and usefulness of chiral perturbation techniques. 
A further step consists on trying to extend this success 
into the strangeness sector. The decay r/ — ► 7r°7r 7r° is 
however not in good agreement with the experimental 
data, see for instance [4|. There can be several sources 
for this discrepancy: i) The f] — rj' mixing may play 
a major role in r\ decays. However it was shown that 
within the chiral framework these effects are small for 
the decay 77 — > 37r ii) Final state interactions are 
significant in this channel Hi) The weight of the 

strange mass, suppressed by large-N c arguments at next- 
to-leading order, can re-emerge via unconstrained higher- 
order counter-terms 0] . Most probably, as shown by the 
unitarization procedure of Q , is some interplay between 
all these issues the cause of the disagreement. Our point 
of view evades a direct calculation in the strong sector 
and is instead based on the experimental data and chiral 
symmetry only. 

Chiral symmetry is explicitly broken by quark masses 
and electric charge. Quark masses are ad hoc parameters 
in the theory and their determination remains one of the 
most challenging tasks. On the experimental side, quark 
masses cannot be measured directly due to confinement. 
One has to relate them to observables which are sensi- 
tive to the variation of their values and try to extract 
the latter from the measurement of these observables. 
In the low-energy regime of the theory, quark masses 
come always multiplied by the quark condensate and one 
can only determine quark mass ratios [lol [ll[ . Limiting 
ourselves to the light-quark sector including strangeness, 
these ratios are equal to one in the case of exact SU (3) 
flavor symmetry. Excluding strangeness, {m u /md) = 1 



in the limit of exact ST/ (2) isospin symmetry. There- 
fore, deviation of the ratio m u /md from unity can be de- 
tected in processes sensitive to isospin symmetry break 
ing. This is the case in the decay 77 — > 37r which is given 
mainly by the difference arid — m u 

The paper is organized as follows. In Section [IT] we 
give the analytic expression of the decay amplitude for 
the process r\ — > 37r° to one-loop order in Chiral Pertur- 
bation Theory (ChPT) including electromagnetic correc- 
tions. In Section [IIII we derive a relation between the up 
and down quark mass difference and the strong part of 
the Dalitz plot parameter. We also give the ChPT predic- 
tion at one loop for the mass difference using data from 
the Dalitz plot. In Section [TVl we use one- loop ChPT cal- 
culation for the decay width and the experimental value 
for the latter to predict the mass difference. Section M 
summarizes our method of extracting the mass difference 
from data only by relating the decay width to the Dalitz 
plot while keeping the strong interaction correction as a 
free parameter. In Section IVII we use our prediction for 
the mass difference to determine the size of the violation 
of Dashen's theorem and up and do wn quark masses. 



II. THE DECAY AMPLITUDE 

Consider the process 

V(p)^7T ( Pl )+7T (p 2 )+7r (p 3 ), (1) 

with Mandelstam invariants 

s=(p- Pl ) 2 , t=(p- P2 ) 2 , u=(p- P3 ) 2 , (2) 
satisfying 

s + t + u = Ml + 3M*o = 3s . (3) 

This is a 3-particle decay. Observables depend on 3 x 
3 free variables. Conservation of 4-momentum reduces 
the number of free variables to 3 x 3 — 4. Since we are 
dealing with bosons, absence of spin implies that in the 
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rest frame of the decaying particle, the orientation of the 
momentum configuration (the angles) is irrelevant and 
hence 3 variables are trivial. There remain 3x3 — 7 
essential variables, s and u, say. 
Let $ define the phase space integral 



$ = /n^ (4) bi+p 2 +P3-p). 

J i=l 1 

The phase space distribution is therefore a constant 

d 2 $ 7T 2 



dsdu AM 2 



(4) 



(5) 



It follows that the density of points in a Dalitz plot is 
proportional to the square of the matrix element. The 
physical region in the (s, u) plane is given by the bound- 
ary of the corresponding Dalitz plot 



AM 2 ^s^ (M v - M w o) 2 , vT < u < 



(6) 



where 



-(3s -s)± ~-a(s, M 2 )X 1/2 (s,M 2 ,M 2 ) , (7) 



and 



a(x, m 2 ) 



1 



4m 2 



(8) 



\(x,m 2 ,ni2) = \x — (m\ — tti 2 ) 2 ] [x — (mi + m 2 ) 2 ](9) 
Integration over u gives 

r 2 r(M n -M„ ) 2 



4M2 , /4M2 



ct(s, M*„ ) A 1/2 (s, M^o , M'*) 



The differential decay rate is defined as 

I 



d 2 T 



2 6 7T 5 M, 



■\M\ 2 d 2 <5>. 



and takes the following form 
d 2 T 1 



|M| 2 , 



(10) 
(11) 

(12) 



by (0. The decay amplitude M. is defined as, 

(ttWIt?) EE i(2lt)H^( V i +p 2 +P3 -p)M , (13) 

and can be deduced from the set of Feynman diagrams 
depicted in Fig. [TJ In order to evaluate these diagrams, 
one has to expand the leading (next-to-leading) order 
chiral Lagrangian to the sixth (fourth) order in the me- 
son field and to first (zeroth) order in the photon field. 
Let q 2 and m q represent generic energy and generic light 
quark mass, respectively. Chiral power counting of the 
diagrams shows that the amplitude counts like q 2 or m q 
at tree level and q 4 , m 2 or m q q 2 at one-loop level in the 
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FIG. 1: Feynman diagrams contributing to the decay ampli- 
tude at one loop. Diagram 1. represents the wavefunction 
renormalization and the -k-t\ mixing. Diagram 2. represents 
the counterterms. Diagram 3. represents the tadpoles. Dia- 
gram 4. represents the unitary corrections. 



isospin limit and in the classic standard scheme, that is, 
q 2 ~ m q . In the isospin breaking case, additional terms 
arise from the two sources of isospin breaking, namely, 
the electric charge and the up and down quark mass dif- 
ference. From the electromagnetic mass of the pion, 



MZ 



0(e 2 ). 



(14) 



we see that e 2 counts as m q . It follows that the additional 
terms count like e 2 or — m u at tree level, and, q 2 (md — 
m u ), m q (m d -m u ), (m d -m u ) 2 , e 2 (m d -m u ), e 2 q 2 , e 2 m q , 
or e 4 at one-loop level. By naive dimensional analysis, 
the e 4 and (md — m u ) 2 terms are suppressed with respect 
to the others and henceforth will be neglected. 

The decay amplitude can be written by symmetry con- 
siderations as 

M(s,t,u) =M(s)+M(t) + M(u). (15) 

The s-channel amplitude can be casted in the following 
general form 

M( S ) = i 1 + M«) + *om(s)} + U'M 

where -Bo is an order parameter for chiral symmetry 
which is related to the quark condensate, 



Be, = — lim 



W) 
o Fl 



(17) 



and F v is the pion decay constant. 

All three corrections in (fTB)) are separately ultraviolet 
finite. Both of them, <J str and J en? , have been already 
considered in the literature, [H, [3] respectively. In ad- 
dition to these we have included a small, with respect 
to the previous ones, term proportional to (m d — m u )e 2 . 
For completeness, we give the analytic expressions for the 
corrections in the Appendix. 
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III. THE DALITZ PLOT 

The Dalitz plot distribution is parameterized as 

\M(s,t,u)\ 2 = \M{s ,s ,so)\ 2 {l + 2az} , (18) 

with so the center of the Dalitz plot and z a dimensionlcss 
parameter defined in terms of pion energies En as 



3Ei - M, 



v 



3 ^ V M„ - iM^o 
and in terms of Mandelstam variables like 



(19) 



1 



{(s-s ) 2 + {t- So ) 2 + (u- So ) 2 } 

w 



2 Af2(M, - 3M.D) 2 

(20) 

The experimental value of a as quoted in the PDG [15[ 



that e — > implies m u — m d . On the other hand, expres- 
sion (f2"3"| relates up and down quark mass difference to 
the electric charge through an observable, namely, a cxp . 

As a first step, we obtain the analytic expressions for 
the a's to first order in the chiral expansion. For com- 
pactness, we use the Gell-Mann-Okubo mass formula, 
AM\ — 3M 2 + M 2 , and obtain the reduced expressions 
for the a's given in the Appendix. 

The main feature of these expressions is the indepen- 
dence on any low-energy constant. This can already be 
seen from the fact that the coefficients of the latter in the 
<5's are polynomials in meson masses and do not depend 
on Mandelstam invariants. Taking the second derivative 
with respect to these invariants gives a null contribution 
to the low-energy constants. The numerical values, at 
"next-to-leading order, of the a's are found to be 



x cxp 



-0.031 ±0.004. 



(21) 



In order to conform to the experimental analysis, we ex- 
pand the square of the amplitude around the center of 
the Dalitz plot with the use of 



j=0 u 



(s - s ) 



d_ 



d 



d 



+(t-s )-^ + (u-s ) — 



du 1 



M(s',t',u')\ (22) 

/ s'—t'—u'—SQ 



We then match with the parametrization (|18p , obtaining 
the scale-independent expression 



B {m d - m u ) 



a CX p — a s tr 



(23) 



Even if we have obtained (|230 at next-to-leading order is 
straightforward to convince oneself that its explicit form 
is maintained at all orders in the strong interaction and at 
first order in the electromagnetic one. Note that the a's 
are given by the second derivative of the corresponding 
<5's, up to normalization factors, evaluated at the center 
of the Dalitz, 

a str = iM2(M„~3M,o) 2 ReC(so), (24) 
a om = ^M 2 (M T ,-3M^) 2 ReS',: Tn (s ), (25) 

= -^M 2 (M V - 3M n o) 2 ReC n ( So ) ■ ( 2Q ) 



Although (|2"3"j) expresses the up and down quark mass 
difference in terms of the electric charge e, it does not 
provide the electromagnetic part of that difference. In 
fact, taking the limit of vanishing e, both numerator and 
denominator in (12311 tend to zero and one cannot claim 



-0.0179, 



(27) 



for the strong piece, and 

a<2 = -0.0011, a<2 = 0.0025, (28) 

for the electromagnetic ones. Notice that (|2"Tj) is roughly 
a factor two bigger and of opposite sign of (|2T)l . Replac- 
ing the a's by their values in (|2"3"1) , we obtain at one-loop 
order, 

B (m d - m„) (1) = -(66 ± 5) MeV 2 . (29) 

As Bo must be positive definite the previous expression 
leads to the conclusion that the up quark is heavier than 
the down quark. The situation does not improve when 
including higher-order corrections or rescattering effects 
in the final state (see Table 5 of [8i]). Moreover a 10% 
correction to the value of (|2"T|) roughly translates in a 4% 
correction to (f2"9")l . A close inspection of (f2"3"|) reveals that, 
disregarding a cm , this inconsistency will always show up 
if a s t r is not bigger, in absolute value, than a eX p- Notice 
that this is against all the findings in the literature tack- 
ling so far the strong sector. We shall attempt in the next 
Sections a different approach using a parametrization of 
available experimental data. 



IV. THE DECAY WIDTH 

As in the preceding section, we can repeat a similar 
procedure for the total decay rate. The decay rate takes 
the following form 

T = B 2 (m d - m u ) 2 (7 tree + 7str + 7cm) 

+B (m d ~m u ) ; y mll , (30) 

where the 7's are obtained, up to a normalization factor, 
by integrating the 5's over the allowed kinematical region, 



7tr 



1 



(31) 



7str = 7^qfReS st r(s) , 

2 

7cm = ^j.FRe<W(s), 



7cm 



T Re 4m (s) 



(32) 
(33) 

(34) The operator T is defined as 
I 



4 

(35) 



2 8 7T 3 M3 ./ 4M 2 



d S a( S ,M^)X 1/2 ( S ,M^,M^)f(s) 



(36) 



The experimental value of the rate as obtained from the 
PDG is 

r cxp = (421 ± 22) eV . (37) 

Using for the amplitude (|16p the Born approximation we 
obtain, 

7troc = 22.31 x 1CT 12 MeV~ 3 , (38) 

while using the next-to-leading expression the contribu- 
tions are split as 

7^ = (35.70 ± 14.14) x 1CT 12 MeV -3 , (39) 

in the strong sector and 

7 « = -(0.02 ± 0.34) x 1CT 12 MeV~ 3 , (40) 
7W = -(2.50 ± 2.44) x 10- 9 MeV _1 , (41) 

in the electromagnetic one. Replacing the 7's by their 
values in (|30p . solving for the quark mass difference and 
discarding the negative root we obtain from r exp the 
next-to-leading order value 

B {m d - m„) (1) = (2717 ± 342) MeV 2 , (42) 

which can be afflicted with sizable and unknown higher 
order contributions. 

Hitherto we have followed two different paths to ob- 
tain Bo(md — m u ), leading to incompatible results: i ) 
First we calculated the 5's in (fl"6|) at one-loop in Chiral 
Perturbation Theory. After we evaluated their second 
derivative at the center of the Dalitz in order to obtain 
the a's (|2"4"|) . We matched the latter with the data coming 
from the Dalitz and obtained for the quark mass differ- 
ence the puzzling value (|29|) . ii) As a second approach we 
integrated the very same expressions for the 5's over the 
physical region in order to obtain the 7 's, (l39"H4"Tj) . We 



matched the latter with the data coming from the de- 
cay width and obtained for the difference the value (I42|) . 
Both results, (|2"9"|) and (|4*2")) are in disagreement. This ob- 
servation is a sufficient motivation to seek for a determi- 
nation of these quantities that is as model independent as 
possible. We shall therefore ignore the theoretical deter- 
minations concerning the strong interaction and provide 
a determination merely based on chiral symmetry and on 
a combination of data from the Dalitz plot and the decay 
rate. 



V. RELATING THE DALITZ TO THE WIDTH 

The decay rate can be obtained by applying two equiv- 
alent methods: i) Either integrating the full amplitude, 
expression (|16[) . ii) Or one can instead integrate the 
parametrization for the amplitude (|18[) . By working out 
explicitly both approaches and matching them, one can 
relate the quark mass difference to the Dalitz parameter 
in the strong sector, a s t r without resorting to any definite 
value for 5 s tr- In doing so, we avoid attempting to dig 
out the how to obtain the piece 5 s t r from first principles 
and we leave it as a free parameter. This is justified by 
comparing the relative size of the perturbative values for 
the (Ts at next-to-leading order and the expectation of 
strong rescattering effects at higher order. 

As a first step we find a relation between Re<5 s t r (so) 
and Bo(m,d — m u ). For that purpose we integrate (|18p 
over the whole Dalitz to obtain the decay rate. Using for 
a and T their experimental values given respectively by 
(|2"T|) and (|3"T)) we determine the distribution at the center 

\M(s ,s ,s )\ 2 = (88.67 ± 4.57) x 10" 3 . (43) 

^From the previous value and the expression of the dis- 
tribution at the center, 



\M{s ,s ,s )f 



) {1 + 2Re<5 str (so) + 2Re<5 em (so) } 



%/3> 2 



B (md - m u ) Re^ cm (so) , (44) 
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derived from (fTCl) , we write Re^ s t r (so) m terms of 
Bo(md — m u ). We shall make use of this relation latter 
in order to pull out Re<5 s tr(so) in favor of Bo(m d — m u ). 

For the last step we start from the definition of 7 str in 
terms of <5 s t r in ([50]) . Expanding the latter around the 
center of the Dalitz and integrating we obtain a relation 
between 7 str , a str and Re5 Btr (s ), 



7str = 27 tree R.e(5 s t r (so) + 



9a s 



3F^M^{M V -3M^) 



Equating now the two expressions for 7 str we obtain a 
second (independent) relation between Re^ s t r (so) an d 
Bo{md — m u ). Solving for Bo(m d — m u ) leads to the 
final expression 

(a stI + a cm )B^(m d -m u ) 2 + a em B (m d -m u )+a cxp = 0, 

(46) 

with 

2 Af 2/ A/ , o M , 2 ^ Re ^cm(s)-L(so)] 
ftcm = g M v (M v - 6M v o ) - 



2F 2 

a cm = --^M 2 {M V - 3M„o) 2 



T{s - so) 2 
J 7 Re S em (s) - (5 C m(so) 



T(s - s ) 



ZF^M^Mr, - 3M^) 2 a exp r cX p 
flcxp ~ 9a exp fi(s - so) 2 + M 2 (M„ - 3M l0 ) 2 H 

At next-to-leading order in the chiral counting 



(47) 



145.374 x 10" 



5 (1) 

"cm 



(48) 

2.905 MeV 2 , (49) 
a cxp = (0.601 ± 0.085) x 10 6 MeV 4 . (50) 

Note that the preceding coefficients are low-energy- 
constant independent. This is due to the fact that, in 
their analytic expression, we integrate the real part of 
the S's subtracted at the center of the Dalitz, that is, 
S(s) — 8(so)- Since the coefficients of the low-energy con- 
stants in the £'s have no kinematical dependence, the 
subtraction allows to get rid of the low-energy constants. 
In other terms, the origin of the errors in the a's, which 
is also the case for the a's, is purely experimental. 

We have now the two relations ([23)1 and (|46|) with the 
two unknowns, Bo(m d — m u ) and a s tr- The first relation 
comes from the Dalitz alone. The second relation is a 
combination of data from both the Dalitz and the width. 
The solution of the system, constraining the down quark 
to be heavier than the up quark, gives the following num- 
bers 

a str = -(0.0305 ± 0.004), (51) 
B (m d -m u ) = (4495 ± 440) MeV 2 . (52) 

If instead we use the recent data from KLOE [16( we 
find a s tr = —(0.0265 ± 0.0100) for the strong interaction 
Dalitz parameter and B (m d -m u ) = (4832±999) MeV 2 
for the quark mass difference. 



VI. SURVEY OF APPLICATIONS 

A. Violation of Dashen's theorem 

We consider the square mass difference Ax-A T , which 
can be split into two pieces 



A T = (A K ) str + (A x - A 7r ) c 



(53) 



where Ap = Mp± — Mp . We have used the fact that 
the pion mass difference is essentially of electromagnetic 
origin 17]. The electromagnetic term in (|53[) vanishes 
at leading chiral order by virtue of Dashen's theorem 
[l8| but is subject to corrections from higher orders. 
These corrections are commonly known as the violation 
of Dashen's theorem. In order to estimate the size of the 
violation we shall introduce the parameter k defined by 



(M 



M 



K°j 



Ml a ) 



(54) 



and rewrite the difference (|53|) at next-to-leading order 
in the following ultraviolet finite form 



(i) 



-B (m d -m u ) {1 + A 8 ti + A Gm } + A c 



The origin of the different terms in (|55[) runs parallel 
to those of ([in]): A str and A em where obtained already 
in 0, U^l , while the sub-leading A cm term is new. Nu- 
merically its contribution turns to be around 8% — 20% 
corrections to the A em value depending on the strong 
low-energy constant set we use. At next-to-leading order 
their numerical values are given by 



A str = 0.277 ±0.340, 
A om = -(0.0014 ±0.0055), 
A em = - (194 ± 1517) MeV 2 



(56) 
(57) 
(58) 



but these are subject to the uncertainty in the low-energy 
constants, as is indeed mainly reflected in (|58[) . Compar- 
ing (|55"j) with ([53"]) and using the estimates ([5TH58"]) we 
obtain, 



= 0.85 ± 1.20. 



(59) 



The central value can be considered moderate in front of 
other estimates [2(| although no solid conclusion can be 
obtained in view of the errors. An alternative, and more 
error free, approach consists on replacing in ([53]) the total 
difference by physical masses, the strong piece by its one- 
loop value and deduce the value of the electromagnetic 
piece and, a fortiori, that of k. We find 



1.4 ± 1.3. 



(60) 



Notice that this value of K contains, beside the violation 
of Dashen's theorem, the effect of higher-order strong 
corrections which can be sizeable. 
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B. Light quark masses 

The parameter Bq is defined by Eq. (fl7|) . We take for 
the pion decay constant and for the quark condensate 
in the chiral limit the unquenched lattice determination 
[H, ^ (2GeV) = (76±3)MeV and (qq) (2GeV) = 
-(242 ± 9) 3 MeV 3 , respectively. This leads to the value 
S (2GeV) = (2454 ± 194) MeV where the main source 
of error comes from Fq. Together with Eq. (|52[) . the 
preceding estimation leads to 

(m d -m„)(2GeV) = (1.8 ± 0.2) MeV . (61) 

Now we use for + m u the lattice determination quoted 
in [H, namely, (m d +m u ){2 GeV) = (7.6±1.6) MeV, and 
obtain 

m„(2GeV) = (2.9 ±0.8) MeV (62) 
m d (2GeV) = (4.7 ±0.8) MeV. (63) 

Both results are, within errors, in agreement with the 
more recent lattice results either using domain wall [2l[ 
[22| , improved Wilson [23[ , twisted mass [24[ or staggered 
fermions (25| . 

VII. CONCLUSION 

In this letter we have used chiral symmetry and data 
to obtain a consistent parametrization of the 77 — > 3tt° 



decay. From the theoretical side we have pointed out 
that the calculation of the strong interaction effects are 
the most probable source of disagreement with data. 
We therefore take this contribution as a free parame- 
ter rather than attempting to evaluate it. The electro- 
magnetic contributions are instead small with respect to 
the strong ones and rather well understood. We incor- 
porate them using Chiral Perturbation Theory. Com- 
bining the Dalitz plot parametrization with data we ob- 
tained a relation between the strong interaction param- 
eter a str and the quantity B^m^ — m,,). Following a 
similar procedure for the width we obtained a second re- 
lation between the same quantities. Solving the system 
for both unknowns we obtained a s t r = —(0.0305 ± 0.004) 
and B (m d - m u ) = (4495 ± 440) MeV 2 . As a first ap- 
plication, we estimated roughly the size of the violation 
of Dashen's theorem (|60[) . As a second application, we 
estimated the values of up and down quark masses (|62p . 



VIII. APPENDIX 



We give here for completeness the expressions for the 
5s and the as appearing in the text 



<5str 



1 1 



240F 2 M 2 K \ 16tt 2 



1 



135s 2 - 15s(18M£ + M 2 + 3M 2 ) + 2M£ (16M£ - 49M 2 + 13Af 2 ) 



-320Af 



K 



M 2 
ivi K 



6(2M| + Ml)Li + (4M 2 - + 11M 2 + 3M 2 )L 5 - 12(2M| + M 2 )L 6 + 48(Af£ - 4M 2 )L 7 - 96M 2 L 8 
][2J (150s - 62M| - 27M 2 - 26M 2 )A(M 2 ) - ( _^ f2 J (2M 2 K + M 2 )A(M 2 ) 



+ 



Ml -Ml 

60M| 
M\ - Ml 



(13M* - 75s)(3M£ + Ml) + 168M£ + 49M%M% - 49M* 



(9s - 22M 2 - 3Af 2 + AM'l)s + M'l{2M 2 K + 7M 2 + 3M^ 



ml) 

B(s,Ml,Ml) 



+ l20MlMlB(s, M 2 , M 2 ) - 240M^M|S(s, M 2 , M 2 ) 



60MI 
Ml -Ml 



(9s - 14M 2 - - 4M 2 - 3M 2 )s 



+Ml(8Ml + M 2 + 3M 2 
2e 2 



B(s,Ml,Ml) +15 



9s - (8M 2 K + Ml + 3M 2 ) 



B'(s,Ml,M 2 K : 



12(K! + K 2 ) - 12(K 7 + K s ) - 6(2K 3 - K 4 ) + 8{K 5 + K e ) - 10(K 9 + K w ) 



Ml 



Ml - Ml 



3{2K 3 - Ki) - 2(K 5 + K 6 



480.F 2 Ml Ml \M 2 K 



16tt 2 Ml - Ml 



135(8M£ - hM'lM'l + M*)s 2 - 5(486M^-M 2 - 57M^M'l + 3M« + 9M*M 2 



(64) 



7 



-AbM^M^M 2 - 96M£ + 72M%M')s - 2M 2 M^(520M| r - 909M;M^ 



+ 13M 2 M 2 - 155M£m£ 



AM 



K 



Ml -Ml 



(300s - <6AM 2 K - 104M ^ - 52M 2 )A(M 2 ) 



AMlM 2 



K 

-M*(8M& - 26M 2 - 6M 2 ) 



j (150, - 92M| + 8M 2 - 26M^)A(M|) + 



60siWl 



18s 2 + 2{AM 2 K - 3M 2 - 22M 2 )s 



B'(s,M 2 ,M 2 ) 



f 30sM 2 
\M*-M* 



9(3M£ - M 2 )s 2 



B\s,M 2 K ,M 2 K ) 



-(36M* + MlM 2 K - M* + 9M|M 2 - 3M 2 M 2 )s + 2M%{SM 2 K + Ml + 3M 2 ) 
-15s 3 M 2 (9s - 8M% - Ml - 3M 2 )B"(s, M 2 K , M|)| , 

2 Mi 3(2K ' - « - <* + *> + + *->} + m 

+15(18M| + Ml + 3M 2 )s - M 2 K {Z2M 2 K - 98M 2 + 26M 2 ) 

+2(-75s + 16M| - 49M 2 + 13M 2 )A(M|-) + 15s 2 (-9s + 8M| + M 2 + 3M 2 )B'(s, M|, M|) J 



- 135s 2 



(65) 



(66) 



a (1) - 

"str — 



a. 



(1) 



1 M 2 (M n -3M v ) 2 „ , 

I^ (M 2 -W)(M 2 + 3M 2 ) X ^ { -54( 47r) - 2 (M 2 - M 2 ) 

-72(M 2 + 3M 2 )S(s , M%, M 2 K ) + 12(M 2 - 5M 2 )(5M 2 - 3M 2 )B'(s , M 2 K , M 2 K ) 

— (M 2 - M 2 )(117M^ - 26M 2 M 2 + 21M*)B"(s , M 2 K , M 2 K ) + 2{M 2 - Ml) 2 (3Ml + M 2 ) 2 B< 3 ' (s , M 2 -, M 2 -) 
+72(Af 2 + 3M 2 )B(s , M 2 , M 2 ) - 24(M 2 + 3M 2 )(M 2 - 2M 2 )B'(s , M 2 , M 2 ) 

+2(M 2 - M 2 )(M 2 + 3M 2 )(9M 2 + 2M 2 )B"(s , M 2 , M 2 ) - 6M 2 (M 2 - M 2 )(M 2 + 3M 2 )B"(s , M 2 ,M 2 ) 

+ 12M 2 (M 2 - M 2 )(M 2 + 3M 2 )B"(s , M 2 , M 2 ) } , (67) 

_1 M 2 (M - 3A^) 2 Rc{ _54( 47 r)- 2 (M4 - 3M 2 M 2 + 18M„ 4 ) 

162F 2 M 2 (M 2 — M 2 )(Ml + 3M 2 ) 2 1 V ' V T ^ 

-36M 2 (11M* - 60M 2 M 2 + 9M*)B'{s , M 2 K ,M 2 K ) 

-12M 2 (7M 2 - 3M 2 )(15M^ - 20M 2 M 2 - 3M$)B"(a , M 2 K ,M 2 K ) 

+Ml(M 2 - Ml)(7Ml - 3M 2 )(3M 2 + M 2 )(27M 2 + M 2 )B^ (s , M 2 K , M 2 K ) 

-2Ml(M 2 - Ml) 2 (3Ml + M 2 fB^( So ,M 2 K , M 2 K ) 

-36(M 2 + 3M 2 ) 2 (2M 2 + 3M 2 )B'(s , M 2 , M 2 ) + 12(M 2 + 3M 2 ) 2 (3M* - 4M 2 M 2 - 3M*)B"{s , M 2 , M 2 ) 
-2(M 2 - Ml){Ml + 3M 2 ) 2 (3M 2 + M 2 )B^ (s , M 2 , M 2 ) } , (68) 



1 M^M^-ZM^) 
hAFl Ml + 3M 2 



Re{27(47r)" 



72M 2 B'(s , M 2 -, Mi) + 3(7M 2 - 3M 2 )(3M 2 + M 2 )B"(s , M 2 K , M 2 K ) 
(M 2 - M 2 )(3M 2 + M 2 ) 2 s( 3 )(s , M 2 -, M 2 -) } , (69) 



expressed in terms of one- and two-point functions. In 
dimensional regularization, these take the following form 



J 



A(m 2 ) 



m 



-2A 



16tt 2 



8 



-2A- Ti(l-ln/i 2 ) ' 



16tt 2 



167T2 7 Q 

r 



(ia;ln[a;mj + (1 — x)m\ — x(l — x)p 2 



(70) 



The symbol /i denotes an ultraviolet scale contained in work was partially supported by European Commission 

RTN Programs MRTN-CT-2004-005104 and MRTN-CT- 



A = 



1 



32tt 2 



4- D 



ln(47r) + 1 - 7e 



(71) 



where 7e is Euler's constant and D the spacetime dimen- 
sion. 
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